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ABSTRACT 

We describe how the coherence function — a Fourier frequency-dependent measure of the linear correlation 
between time series measured simultaneously in two energy channels — can be used in conjunction with energy 
spectra, power spectra, and time delays between energy channels to constrain models of the spectrum and vari- 
ability of X-ray binaries. Here we present a procedure for estimating the coherence function in the presence of 
counting noise. We apply this method to the black hole candidates Cyg X-l and GX 339-4, and find that the near 
perfect coherence between low- and high-energy X-ray photons rules out a wide range of models that postulate 
spatially extended fluctuating emission, thermal flares, and overlapping shot-noise. 

Subject headings: black hole physics — methods: statistical — X-rays: stars 



> 

in 

(N 

o 

so : 
os ■ 

Or 
i 1 

o ■ 
a: 



1. INTRODUCTION 

Rapid (T < 1000 s) aperiodic variability is common to all 
types of X-ray binaries. It arises in the immediate vicinity of 
the compact object and provides a probe of changes in physi- 
cal parameters such as accretion rate, optical depth, and tem- 
perature. Power spectra and cross spectra have been applied 
widely to study variability in X-ray binaries as a function of 
frequency and energy, as well as to measure time delays be- 
tween intensity variations at different energies (cf. van der Klis 
1989). Here we describe the coherence function — a measure 
of the degree of linear correlation between two time series as 
a function of Fourier frequency — and how this provides strong 
constraints on models of X-ray binaries, especially recent theo- 
retical models that attempt to correlate energy spectra with ape- 
riodic variability behavior (cf. Nowak 1994, Miyamoto et al. 

1994, Miller 1995, Nowak & Vaughan 1996). We expect the 
coherence function to be widely applicable to Rossi X-ray Tim- 
ing Explorer (RXTE) observations. 

No current model mimics all facets of a system as complex 
as an X-ray binary. However, spectral and dynamical models 
should qualitatively reproduce — or at the very least not funda- 
mentally disagree with — all of their most basic average prop- 
erties: the energy spectrum, power spectrum (PSD), and time 
delays and coherence between different energies. (For ways in 
which a spectral model can disagree with timing data, cf. Miller 

1995, Nowak & Vaughan 1996.) Most theories model the en- 
ergy spectra, fewer model the PSD, fewer still consider phase 
lags, and, to date, none consider the coherence function. Here 
we describe several examples that lead to unity coherence, as 
well as several that generically lead to a loss of coherence. 

We apply the coherence function to the black hole candidate 
Cyg X-l in its low (hard) state and GX 339-4 in its very high 
state. The former shows unity coherence over a wide range of 
Fourier frequencies and energy bands, whereas the latter shows 
a sharp drop in coherence between low- and high-energy bands. 
We present these data to show the application of the statistical 
methods to real data, and to show that there are intrinsic differ- 
ences among distinct physical states and systems. 



2. COMPUTING THE COHERENCE FUNCTION 

Let xi(k) and x 2 (k) be the X-ray flux observed simultane- 
ously in two energy channels at times tk. The Fourier trans- 
forms ofx\{k) and x 2 (k) at frequency Vj areXi(j) andA^j), 
respectively. We will usually use lowercase variables to in- 
dicate time series and uppercase variables to indicate Fourier 
transforms. From the Fourier transforms one can construct the 
power spectra [P x (j) = |X x (j)| 2 and P 2 (j) = \X 2 (j)\ 2 ] and 
the (complex-valued) cross spectrum [C(j) = Xl(j)X2(j)], 
whose argument is the phase difference between intensity fluc- 
tuations in the two channels at frequency Vj which can be con- 
verted to a time delay St(j) = arg[C(j)]/ (2m>j). 

Measured X-ray light curves always contain noise. For most 
X-ray observations, the dominant noise source is counting, or 
Poisson, noise. We denote signal by s and noise by n, and write 
xi(k) = si(k) +ni(k), Xi(j) = Si(j)+Ni(J), and likewise 
for channel 2. Hereafter we usually drop the explicit frequency 
dependence to simplify notation. For power spectra, Poisson 
noise adds to the signal, and Pi = \Si\ 2 + |Ai| 2 . The signal 
power, | Si | 2 can be estimated by subtracting the Poisson noise, 
| TVi | 2 , from the measured power. 

We now consider the idealized case of two simultaneous ran- 
dom processes, 1 and 2, and let Si(i) and S2(t) be noiseless 
signals drawn from these processes. It is usually possible to 
find a linear transform, h(r), called the transfer function, that 
relates Si(t) and S2(t) via 



«a(t) 



h(t — t)s\(t) dr, 



(1) 



or, equivalently, 82(f) = H(f)Si(f), where S and H are 
the Fourier transforms of s and h. If H(f) is the same for 
all realizations of the two processes, the processes are said 
to be coherent at frequency /. In that case, | (C(/)) | 2 = 

l^(/)| 2 <|5i(/)| 2 ) 2 = (|Si(/)| 2 )(|S 2 (/)| 2 >, where angle 
brackets denote an average over an infinite set of signals from 
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the same processes. The coherence function is then defined as 

I <C(/)> I 2 



(|^i(/)| 2 )(|5 2 (/)| 2 )' 



(2) 



where we affix a subscript / when computing the intrinsic co- 
herence between noiseless signals. It is clear that when Si(f) 
and S 2 (/) are related by a linear transform, vf (/) = 1. Co- 
herence is a measure of the fraction of the rms amplitude of 
one process at a given frequency that can be predicted from the 
other by a linear transform (Bendat & Piersol 1986). It is im- 
portant to note that unlike powers and time delays, coherence 
can only be computed for an ensemble of independent mea- 
surements. For noiseless measurements from a finite number, 
m, of independent samples, the statistical uncertainty of the co- 
herence function is £7/ = (2/to) 1 / 2 (1 — T/)/|Tz| (Bendat & 
Piersol 1986). 

The coherence function of real data is given by analogy to 
equation (^|), except that angle brackets are used to denote 
an average over a finite number, m, of independent measure- 
ments. The coherence function of noisy data will always be 
less than unity. The intrinsic coherence can be estimated by 
correcting each term in the discrete analogy to equation (fy for 
counting noise. For the powers in the denominator, one sets 
|Si| 2 = P x - liVxl 2 , and similarly for \S 2 \, where P x and P 2 
are the measured, noisy power spectra of the signals. 

Writing out | (C) | 2 in terms of signal and noise yields 



in which case (a 2 ) = (s 2 ) + (n 2 ) and Sa 2 = \[2r? . The 
intrinsic coherence can be usefully estimated when the follow- 
ing conditions are met: s > n, |S'i| 2 /|iVi| 2 > 1/y/rri, and 

|S 2 | 2 /|iV 2 | 2 > i/y^. 

The following is a recipe for estimating the value and un- 
certainty of the intrinsic coherence from measurements of (C) , 
Pi and P 2 . Average powers and cross spectra should be con- 
structed from unnormalized measurements. In practice it will 
often be necessary to correct the Poisson noise level for instru- 
mental dead time (cf. van der Klis 1989 for a detailed discus- 
sion). Dead time can introduce correlations between energy 
channels that enhance or diminish coherence. These are intru- 
ment dependent and beyond the scope of this Letter. The terms 
high power and high measured coherence used below denote 
powers that satisy l^ 2 greater than a few times |iV| 2 / \pm in 
each channel and measured coherence that satisfies greater 
than a few times t? j {P\P 2 \ respectively. In many cases, these 
conditions will be satisfied at some frequencies but not at oth- 
ers. 



High powers, high measured coherence.- 



7? 



(C) 



\Si\ 2 \S 2 \ 2 



1 ± m 



-1/2 



2n m 



(\{C) 



\Si\ 



mSj 2 



1/2 



(8) 



(C) 



(StS 2 ) + (SIN 2 ) + (N*S 2 ) + (N*N 2 ) | 2 . (3) 



The first term is the average cross spectrum of the signals. Pois- 
son noise in one channel is uncorrected with Poisson noise 
in the other channel and with signal in either channel. Thus, 
I (C) 1 2 is the squared magnitude of the sum of SIS 2 and a 
complex random variable, <r, whose real and imaginary parts 
each have zero mean and variance 



var(^) 



1 



\Si 



\ 2 \N 2 \ 2 - 



\Ni 



\ 2 \S 2 \ 2 - 



We define s 2 

(\Si\ 2 \N 2 \ 



|iVi| 2 |iV 2 | 2 )/m . (4) 

(SIS 2 ) | 2 , a 2 = I (C) | 2 , and n 2 = 
K| 2 |/V 2 | 2 )/m. The probability 



distribution of the measured cross spectrum, | (C) | 2 , given an 
intrinsic cross spectrum and Poisson noise, is equivalent to the 
distribution of total power in an individual frequency bin of a 
noisy power spectrum; it is given by 



P(a 2 \s 2 ,n 2 ) 



n 2 exp 



-(a 2 



2as 
n 2 



(5) 



This is the Gaussian limit, and in essence it corresponds to opti- 
mally filtering the measured coherence. The first three terms in 
the uncertainty come from uncertainties in Poisson noise, and 
the last term is from the statistical uncertainty in the intrinsic 
coherence. 

High powers, low measured coherence. — Use equation (||) 
to determine confidence limits, | (S1S 2 ) ^ and | (SiS 2 ) | 2 nax , 
on I (SIS 2 ) \ 2 . Confidence limits on j 2 are then 



(sts 2 



\Si\ 2 \S 2 \ 2 



1 ± 



|JVi 



m\Si\ 



\S 2 



if 



(9) 

Low powers. — This difficult case arises in weak sources and 
at high frequency. In practice the la errors in this case are 
likely to extend nearly from to 1 . The authors do not know a 
closed form for p(7 2 ). We recommend using equation (||) and 
the Gaussian probability distributions of |Si| 2 and \S 2 \ 2 to em- 
perically map out p(j 2 ) and determine confidence limits on the 
intrinsic coherence. 



(Groth 1975, Goodman 1985), where Iq is the zero order Bessel 
function of the first kind. This can be inverted to yield the prob- 
ability distribution of s 2 given a measured cross spectrum and 
Poisson noise, 



I 2i 2 2\ 

p{s \a ,n ) 



(ns)~ 



■ exp 



2s 2 ) 



2n 2 



I (2as/n 2 ) 
I {a 2 /2n 2 )' 



(6) 

(Chakrabarty 1996). For a ^ n, equation || can be approxi- 
mated by the Gaussian 



P(a 2 \s 2 ,n 2 ) 



2nay / 7r 



exp 



-(a- 



(7) 



3. EXAMPLES OF INCOHERENT SOURCES 

Thermal Flares. — Let us imagine that the observed variability 
is completely the result of local temperature fluctuations with 
some time dependence, T(t), which can be arbitrarily compli- 
cated so long as it is statistically stationary. Let the intrinsic 
time series, si(t) and s 2 (t), be photon count rates in narrow 
frequency bands (yi, v\ + 5vi) and {v 2 , v 2 + Sv 2 ). The ob- 
served count rate in frequency band 1 is approximately pro- 
portional to v 2 5i>i (exp [hvi/kT{t)\ — 1)~ , and likewise for 
energy band 2. If both hv\ and hv 2 <C kT(t), then we are 
on the Rayleigh- Jeans portion of the spectrum, in which case 
s\(t) (x s 2 {t) oc T(t). Thus there is a linear transfer function 
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(a constant) between the two channels, and therefore the coher- 
ence function will be unity. On the other hand, if kT(t) <C hv\ 
and hv2, then we are on the Wien tail of the spectrum, in which 
case we have s 2 {t) oc s 1 (i)( </2 / 1 ' 1 ) oc exp[-hv 2 /kT(t)]. Thus 
there is a nonlinear transfer function. In general, such a transfer 
function will take power from a frequency / and distribute it 
among harmonics of / (cf. Bendat & Piersol 1986) and there- 
fore lead to a loss of coherence. 

Multiple Flaring Regions. — A number of models for X-ray 
variability in Galactic black hole candidates associate different 
timescales with different physical regions of an accretion disk 
(cf. Nowak 1994, Miyamoto et al. 1994). If more than one 
region contributes to the signal in both energy bands, then it is 
possible for the coherence function to be less than unity, even 
if individual regions produce perfectly coherent variability. 

Consider two flaring regions: one produces a time series 
q(t) = qi(t) + g 2 (t), and the other produces a time series 
r(t) = 7*i (t) + r2(t) — each in energy bands 1 and 2, respec- 
tively. Furthermore, assume that there is a constant linear trans- 
fer function that relates qi(t) to q2 (t), as well as another con- 
stant, linear transfer function that relates r±(t) to r 2 (i), with 
qi (t) and r% (t) otherwise being completely uncorrected. De- 
noting the intrinsic time series observed in band 1 as Si(i) = 
qi (t) + r\ (i), and the intrinsic time series observed in band 2 as 
S2 (t) — q2 it) + T2 {€), the intrinsic coherence function between 
the two bands is then 



7? 



Q?Ql + J Rf J Ri + 2|Qi||Q 2 ||i?i|| J R 2 |cos(<J0 r 



Q1Q2 + ^1^2 + Ql^-2 + Q2R1 



where 59 q and 59. 



(10) 

are the mean Fourier phase differences be- 
tween Qi(f) and Q2(f) and between i?i(/) and R 2 (f), re- 
spectively. (The quantities on the right hand side of equation 
( |l0| ) refer to their mean values.) One can easily show that for 
this case < 1, and is only equal to unity if both 59 q = 69 r 
and IQ1I/IQ2I = l-Ril/l-Ral- That is to say, the coherence func- 
tion is unity if and only if the same linear transfer function that 
takes qi(t) to q2(t) also takes to r 2 (i). 

4. EXAMPLES OF COHERENT SOURCES 

Here we present an idealized model that produces unity co- 
herence while producing phase lags qualitatively similar to 
those seen in Cyg X-l (i.e. phase lags approximately inde- 
pendent of Fourier frequency; Miyamoto et al. 1992). Imag- 
ine that we have a disk with a stochastic source of (linear) 
surface density perturbations at the center that propagate out- 
wards at some constant speed c s . If we take the source to be 
given by (Aitr) -1 5(r)a s (t), then the resultant surface density 
perturbations — <r p (r,t), at radius r and time t (using the two- 
dimensional wave propagation Green's function; cf. Morse & 
Feschbach 1953) — are given by 



a p (r,t) = / dt'a s (t') 



(t-f) - r/c s 



7c? 



(ii) 



where 9 is the step function. Let the fluctuating signal in band 
1, si(t), be equal to the surface density perturbation times 
a weighting function, gi{r), integrated over the entire disk. 
The radially integrated Green's function is just a linear trans- 
fer function, thus the Fourier transform of si (t) becomes 

Si(f) = E s (/) J *2n 2 r 9l (r) (^^) dr , (12) 



which we define as E s (/)Ti (/), and where Hq (x) is the Han- 
kel function of x as well as the Fourier transform of the two- 
dimensional Green's function. Similarly we have S%{f) = 
E S (/)T 2 (/). As both S\ and S 2 are linearly related to the 
source E s , there is a linear transfer function equal to T 2 /Ti 
that takes S\ to 5 2 , and therefore there is unity coherence. We 
could write such a transfer function only because the driving 
source was separable in time and space (in this case being spa- 
tially localized). 

Equation ( |12| ) is analytically tractable if we set gi(r) — 
exp(— atir/ro) and 52 (^) = exp(— a^r/ro), as one might have 
for weighting functions representing flux from a Wien tail. As 
will be discussed elsewhere (Vaughan et al. 1997), the result- 
ing time lags are qualitatively similar to those seen in Cyg X-l. 
We require ro/c s ~ 1 to obtain quantitative agreement with 
the observations, which implies an extremely slow propagation 
speed, c s ~ 10~ 4 c. This mechanism is qualitatively similar to 
that proposed by Manmoto et al. (1996), where thermal waves 
are launched from large disk radii toward the center, and then 
reflect from the disk inner edge as acoustic waves. However, 
the mechanism of Manmoto et al. (1996) cannot preserve co- 
herence if, as in the example of §3, there are multiple, spatially 
distributed source functions (in the X-ray -emitting regions) for 
the waves. Furthermore, the acoustic waves shock upon reflec- 
tion, which is an inherently non-linear process that also will 
lead to a loss of coherence, if emission from the reflected shock 
is an appreciable fraction of the total observed emission. 

Adding additional mechanisms in the form of successive 
transfer functions will lead to unity coherence so long as each 
component is itself linear. In the above example, we could have 
invoked Comptonization as a subsequent transfer mechanism, 
which would have introduced additional time delays. So long 
as the Compton cloud is static, coherence will be preserved (cf. 
Nowak & Vaughan 1996). 



5. OBSERVATIONAL EXAMPLES 

Of the above physical situations where coherence is pre- 
served and those where it is not, which cases are observed in na- 
ture? Surprisingly, the answer seems to be that unity coherence 
is the norm, despite the variety of ways of weakening coher- 
ence. Here we present two black hole candidates: GX339-4 in 
its very high state, and Cyg X-l in its low state. Both have been 
modeled with disks plus Compton coronae (cf. Miyamoto et al. 
1991; Dove, Wilms, & Begelman 1996 and references therein), 
although there are differences between them: the very high state 
has a substantial soft (~ 1 keV) component, the low state has 
none; the hard tail of the very high state is softer than the hard 
tail of the low state (~ 2.5 photon index compared to ~ 1.7); 
the very high state has lower variability (rms ~ 10 — 20%) than 
the low state (rms ~ 40%). Yet, as shown in Figure 1, both sys- 
tems apparently have unity intrinsic coherence, at least between 
low energies and at low frequency. 

Figure 1 presents coherence plots derived from Ginga data. 
The methods described above have been used to filter the noise 
from the data. Poisson noise and uncertainties in the Ginga 
dead time make the coherence estimates unreliable above ~ 10 
Hz ; however, at frequencies below 10 Hz the intrinsic coher- 
ence is essentially unity. This result holds true for Cyg X-l over 
a wide range of energy channels. On the other hand, the GX 
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339—4 variability (the PSD is shown as Fig. 4b in Miyamoto 
et al. 1991) shows unity coherence between low energies but 
shows a dramatic drop in coherence between the 2.3 — 4.6 keV 
and 13.9 - 37.1 keV channels. 

Considering the mechanisms discussed above, this loss of co- 
herence is not surprising. Perhaps the 13.9 — 37.1 keV band 
corresponds to a Wien tale, and, therefore, to a nonlinear trans- 
fer function. It is more difficult to explain the observed unity 
coherence. We expect that most shot models (cf. Lochner et 
al. 1991) will produce less than unity coherence. For exam- 
ple, the kinematic model of Nowak (1994) for the very high 
state of GX 339^- successfully reproduced the observed PSDs 
and phase lags. However, as the viscous and thermal fluctua- 
tions were distributed over a large range of radii and overlapped 
in time, the mechanism led to less than unity coherence (not 
shown in that paper) between all energy bands. In order to pre- 
serve coherence, we always require linear responses, and we 
usually require the following: localized sources and/or local- 
ized responses; a (temporally) uniform source throughout the 
disk; or a uniform response throughout the disk. Most of these 
features are absent from current models. 

The data described above will be presented in greater de- 
tail in Vaughan et al. (1997). Here we have presented them 
to demonstrate that there are observational and, by extension, 
physical differences between the systems, despite the fact that 
they are both commonly fit with coronae models. The coher- 
ence function therefore offers us another tool to help distinguish 
among theoretical models. 

6. CONCLUSIONS 

We have described a statistic, the coherence function, that 
is derivable from the cross spectrum used to compute phase or 



time lags between two times series. This statistic has tradition- 
ally been ignored; however, it contains additional information 
about the system, and the methods described above can mini- 
mize the noise effects. Whenever one calculates the phase lag, 
one can and should also calculate the coherence function. 

In general, coherence will be lost whenever there is a nonlin- 
ear transfer function between two channels, or whenever there 
are multiple, uncorrected (linear or nonlinear) transfer func- 
tions between two channels. However, our experience has been 
that coherence is preserved more often than not. This is a chal- 
lenge for theoretical models, as there are many more mecha- 
nisms for destroying coherence than there are for preserving 
coherence. 

We presented data from two black hole candidates: GX 339- 
4 in its very high state and Cyg X-l in its low state. At low fre- 
quency and low energy both preserve coherence, although GX 
339-4 loses coherence between high- and low-energy bands. 
Theoretical models have addressed observations of phase and 
time lags (cf. Nowak 1994, Manmoto et al. 1996, etc.), as 
measured with EXOSAT, Ginga, and other X-ray timing instru- 
ments. No model currently accounts for the coherence proper- 
ties detected in GX 339-4 and Cyg X-l. Forthcoming RXTE 
observations, however, coupled with the tools presented above, 
can help provide new insights into the physical mechanisms at 
work in these systems, as well as new challenges for theorists 
to meet. 

The authors would like to acknowledge useful conversa- 
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port from NASA grants NAG 5-3239 (B.A.V.), NAG 5-3225 
(M.A.N.), and the Netherlands Organization for Scientific Re- 
search (NWO) under grant PGS 78-277. 



REFERENCES 



Bendat, J & Piersol, A. 1986, "Random Data; Analysis and Measurment Pro- 
cedures", (New York: Wiley and Sons) 
Chakrabarty, D. 1996, PhD Thesis, Caltech 
Dove, J., Wilms, J., & Begelman, M.C. 1996, ApJ, submitted 
Goodman, J.W. 1985, "Statistical Optics" (New York: Wiley) 
Groth, E.J., 1975, ApJS, 29, 285 
Lochner, J.C. et al. 1991, ApJ, 376, 295 

Manmoto, T., Takeuchi, M., Mineshige, S., Matsumoto, R., & Negoro, H. 1996, 

ApJ, 464, L135 
Miller, M.C. 1995, ApJ, 441, 770 

Miyamoto, S., Kimura, K., Kitamoto, S., Dotani, T., & Ebisawa, K. 1991, ApJ, 



383, 784 

Miyamoto, S., Kitamoto, S., Iga, S., Negoro, H., & Terada, K., ApJ 1992, 391, 
L21 

Miyamoto, S., Kitamoto, S., Iga, S, & Hayashida, K. 1994, ApJ, 435, 398 
Morse, P.M. & Feschbach, H. 1953, Methods of Theoretical Physics (New 

York: McGraw Hill) 
Nowak, MA. 1994, ApJ, 422, 688 

Nowak, MA. & Vaughan, BA. 1996, MNRAS, 280, 227 

van der Klis, M. 1989, in "Timing Neutron Stars", ed. H. Ogelman and E.P.J. 

van den Heuvel (Dordrecht: Kluware Academic Publishers), p. 27 
Vaughan, BA. et al. 1996, in preparation 



5 



Vaughan & Nowak 



Vol. 999 



o 
o 



0.01 0.10 1.00 10.00 



2 - 



0.8 - 




OA - 







GX339-4 



















0.01 0.10 1.00 10.00 
Frequency (Hz) 



- 1.2 
1 .0 

- 0.8 
0.6 

- 0.4 

0.2 
0.0 



Jl 



0.01 0.10 1.00 10.00 



Fig. 1. — Coherence plots derived from Ginga data of Cyg X-l in its low state (left: 1 — 2.3 vs. 9.2 — 11.5 keV), and GX 339^1 in its very high state (middle: 
2.3 — 4.6 vs. 4.6 — 9.2 keV; right: 2.3 — 4.6 vs. 13.7 — 37.1 keV). The high-frequency rollover in Cyg X-l is likely due to a misestimation of the Poisson noise level. 
Coherence exceeding unity in GX 339^t may be the result of dead time effects. The rapid coherence rise on the right is at the quasi-periodic oscillation frequency 
(cf. Miyamoto et al. 1991). 



